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PREFACE 


Under  certain  conditions  of  flow  near  a  slender  body,  when  surface 
injection  rates  are  large  and  are  not  dependent  on  convective  heating, 
it  may  be  possible  to  "blow  off"  the  boundary  layer.  The  flow  field  of 
such  a  body  can  then  be  investigated  using  an  inviscid  flow  model,  if 
the  resulting  shear  layer  is  sufficiently  thin. 

This  Memorandum  considers  a  special  class  of  such  problems  called 
"Blowhard  Problems,"  where  the  injectant  layer  is  thicker  than  the  shear 
layer,  but  thin  enough  to  preserve  the  slenderness  of  the  effective 
body  produced  by  the  injection.  High  mass- transfer  rates  might  be 
achieved  by  transpiration  through  porous  walls  or  by  ablation  due  to 
radiation  heating.  The  results  of  this  Memorandum  should  be  useful  in 
the  interpretation  of  wind  tunnel  tests  and  in  the  prediction  of  aerody¬ 
namic  or  re-entry  performance  under  conditions  of  high  mass  transfer. 

One  of  the  authors,  J.  D.  Cole,  is  a  professor  at  the  California 
Institute  of  Technology,  and  a  consultant  to  The  RAND  Corporation. 


ABSTRACT 


A  simple  model  is  presented  for  the  analysis  of  inviscid  flow 
fields  over  slender  bodies  accompanied  by  high  rates  of  surface  mass 
transfer.  This  model  assumes  a  thin,  inviscid  injectant  layer,  which 
is  separated  from  the  outer  flow  by  a  contact  discontinuity.  Inviscid 
boundary-layer  equations  are  shown  to  be  applicable  within  the  injectant 
layer,  and  the  pressure  gradient  is  found  by  matching  pressure  and  flow 
direction  at  the  dividing  streamline.  Using  stream-function  variables, 
an  Abel  integral  equation  is  derived,  and  the  inverse  problem  of  de¬ 
termining  the  injection  distribution  from  the  shape  of  the  dividing 
streamline  is  solved  analytically  for  supersonic  and  hypersonic  flow 
over  flat  plates  and  cones. 
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I .  introduction'' 


The  problem  considered  in  this  Memorandum  is  the  construction  of  a 
model  of  flow  past  bodies  with  high  rates  of  mass  addition  or  hard  blow¬ 
ing  at  the  surface.  Blowing  is  considered  "hard"  when  the  velocity  com¬ 
ponent  normal  to  the  wall  is  an  order  of  magnitude  larger  than  that 
usually  expected  in  a  boundary  layer  with  no  blowing.  Under  such  cir¬ 
cumstances  the  boundary  layer  can  be  expected  to  blow  off  the  body  and 
become  a  free  shear  layer  separating  the  blown  flow  from  the  free- stream 
flow.  As  a  limiting  case  corresponding  to  high  Reynolds  numbers,  this 
free  shear  layer  can  be  regarded  as  a  slipstream,  and  the  inviscid  flow 
outside  the  slipstream  can  be  matched  to  the  inviscid  flow  within  it. 

The  flow  between  the  slipstream  and  the  wall  is  essentially  rotational. 
Under  various  simplifying  assumptions,  solutions  can  be  obtained  for 
this  inner  flow  layer  and  matched  to  suitable  external  flows.  In  this 
section  we  present  some  details  of  the  general  concepts  employed.  In 
Section  II  the  basic  approximation  and  the  solution  are  given  for  the 
simplest  case.  Additional  special  cases  appear  in  Sections  III  and  IV, 
and  the  general  conclusions  are  presented  in  Section  V. 

It  is  well  known  that  the  structure  of  both  turbulent  and  laminar 
boundary  layers  can  be  changed  dramatically  by  fluid  injection  from  the 
boundary.  Injection  reduces  skin  friction  and  heat  transfer,  increases 
boundary -layer  thickness,  and  destabilizes  the  laminar  boundary  layer. 
Theoretical  evidence ^  indicates  that  the  laminar  boundary  layer  on  a 
flat  plate  "blows  off"  the  wall  when  the  blowing  rate  is  such  that 


P  v 
w  w 

P  u 

CO  CO 


Injection  rates  higher  than  this  cannot  be  treated  by  boundary- layer 
theory. 


Some  highlights  of  this  Memorandum  were  presented  orally  at  the 
Fluid  Dynamics  Division  of  the  American  Physical  Society  Meeting  in 
Hawaii,  September  1965. 
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(2) 

Some  experimental  evidence'  indicates  that  the  low- speed  turbu¬ 
lent  boundary  layer  on  a  flat  plate  "blows  off"  at  an  injection  rate 
Pwvw/PmUra  approximately  .02.  Within  classical  flat-plate  laminar 
boundary- layer  theory,  blow-off  appears  as  a  singularity;  it  is  not 
possible  to  obtain  solutions  to  the  boundary- layer  equations  that  sat¬ 
isfy  the  given  boundary  conditions  beyond  the  blow-off  point.  Reference 
3  contains  a  discussion  of  the  effects  of  mass  injection  on  incompres¬ 
sible  boundary  layers  and  an  investigation  of  the  blow-off  singularity 
on  a  flat-plate  boundary  layer  with  a  constant  injection  rate, 

(4) 

Pretsch  examined  the  Falkner-Skan  flows  for  asymptotically  high 
injection  rates  and  favorable  pressure  gradients.  He  demonstrated  that 
solutions  are  possible  for  arbitrarily  large  injection  rates,  and  showed 
that  skin  friction  approaches  zero  as  j3/(-f  ),  where  (f^)  is  the  wall 
injection  rate  in  the  appropriate  similarity  variable  and  3  is  the  pres¬ 
sure-gradient  parameter.  Pretsch  also  showed  that  neglecting  the  vis¬ 
cous  terms  in  the  equation  of  motion  for  high  injection  rates  produces 
an  inviscid  balance  between  inertia  and  pressure  gradient.  The  result¬ 
ing  inviscid  equation,  solved  by  Pretsch,  has  a  discontinuity  in  vortic- 
ity  at  the  dividing  streamline.  In  addition,  for  certain  values  of  3, 
the  shear  becomes  infinite  at  the  dividing  streamline,  making  it  neces¬ 
sary  to  introduce  a  viscous  shear  layer  at  the  dividing  streamline  to 
resolve  this  singularity.  An  extension  of  Pretsch' s  original  inviscid 
analysis  to  the  details  of  the  shear  layer  is  given  in  Refs.  5  and  6. 

More  recently,  Libby^^  and  Vinokur ^  have  discussed  the  problem 
of  the  stagnation  region  of  a  sphere  in  hypersonic  flow  when  mass  is 
injected  from  the  surface  of  the  sphere.  Their  papers  treat  the  prob¬ 
lem  of  inviscid  flow  with  mass  addition  within  the  framework  of 
hypersonic,  constant-density,  blunt-body  flow.  Their  solutions  are 
characterized  by  two  regions  separated  by  a  surface  of  discontinuity  -- 
the  dividing  streamline.  The  inner  region  consists  of  mass  injected 
from  the  wall,  and  the  outer  region  of  fluid  that  has  passed  through 
the  bow  shock  wave.  The  inner  flow  may  be  rotational  or  irrotational , 
depending  on  the  manner  of  fluid  injection,  whereas  the  outer  flow  is 
rotational  because  of  the  curvature  of  the  shock  wave.  Since  there  is 
a  stagnation  point  in  the  flow  through  which  the  dividing  streamline 
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passes,  the  total  pressure,  static  pressure,  and  flow  direction  are 
continuous  as  the  discontinuity  is  approached  from  either  the  inner 
injection  zone  or  the  outer  shock  layer.  Vinokur  found  that  the  thick¬ 
nesses  of  the  shock  layer  and  the  injection  layer  are  functions  of  the 

2  2 

nondimens ional  injection  momentum  p  V  / (p  V  )  alone. 

w  W  co  co' 

Ting^  '  has  recently  considered  a  model  of  distributed  injection 
on  a  flat  plate,  with  the  restriction  that  the  inner  flow  is  incompres¬ 
sible  and  potential  and  the  outer  flow  is  linearized. 

In  this  Memorandum  we  consider  the  effects  of  mass  transfer  on 


flows  that,  in  the  absence  of  mass  transfer  or  viscous  effects,  would 
have  a  pressure  gradient  of  zero  along  the  surface.  Examples  of  such 
flows  are  supersonic  wedge  and  cone  flows.  For  such  flows,  high  injec¬ 
tion  rates  could  result  in  blow-off,  making  it  necessary  to  re-evaluate 
the  usual  procedures  of  boundary  layer  theory. 

According  to  boundary- layer  theory  (as  some  typical  Reynolds  num¬ 
ber  approaches  infinity),  the  inviscid  flow  around  the  object  is  first 
calculated,  establishing  a  velocity  and  a  pressure  field  for  the  bound¬ 
ary  layer.  This  inviscid  flow  typically  slips  along  the  solid  surface, 
and  a  viscous  boundary  layer  along  the  surface  must  be  introduced  to 
satisfy  a  no- slip  condition.  In  this  theory,  blowing  must  vanish  as 
Re  ->  00 .  Higher-order  interaction  effects,  such  as  flow  due  to  displace¬ 
ment  thickness,  can  then  be  calculated.  When  the  blowing  rate  is  suf¬ 
ficiently  high  this  approximation  scheme  may  not  be  suitable.  The  first 
step  in  a  more  appropriate  calculation  scheme  is  to  find  the  inviscid 
flow  past  the  object  with  blowing.  This  flow  contains  a  dividing 
streamline,  off  the  body  surface,  which  separates  the  injected  fluid 
from  the  oncoming  stream  (see  Fig.  1).  In  inviscid  theory,  there  is 
a  jump  in  tangential  velocity  across  the  separating  streamline.  The 
flow  in  the  injection  layer  is  generally  rotational,  but  may  to  the 
first  approximation  be  considered  inviscid.  The  viscous  boundary  layer 
becomes  a  free  shear  layer  which  resolves  the  discontinuity  at  the 
dividing  streamline. 

In  the  following  sections  we  consider  a  simple  case  of  this  type 
of  flow,  in  which  blowing  produces  only  a  thin  region.  An  assessment 
of  orders  of  magnitude  shows  that  the  flow  in  the  blown  layer  can  be 
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described  by  inviscid  boundary-layer  equations.  These  equations  are 
simple  enough  to  provide  an  explicit  relation  between  the  distribution 
of  blowing  velocity,  the  distribution  of  pressure,  and  the  shape  of 
the  dividing  streamline.  If,  in  addition,  there  is  some  simple  rela¬ 
tion  between  shape  and  pressure,  as  for  example  in  linearized  super¬ 
sonic  flow,  then  explicit  values  of  all  the  quantities  can  be  found. 

This  simple  theory  should  be  applicable  whether  the  shear  layer 
is  laminar  or  turbulent,  so  long  as  it  is  thin  relative  to  the  region 
of  blown  flow. 

It  should  also  be  noted  that  the  case  of  strong  suction  is  com¬ 
pletely  different  from  strong  blowing;  the  boundary  layer  always  oc¬ 
curs  at  the  wall.  However,  since  the  boundary  layer  does  not  occur  on 
a  streamline  of  the  inviscid  flow,  the  boundary-layer  equations  are 
ordinary  differential  equations  and  the  boundary-layer  thickness  is 
0 (1/Re) . 


II.  INCOMPRESSIBLE  INNER  LAYER 


Consider  flow  past  a  flate  plate  (x  >0)  in  a  uniform  stream  with 
velocity  U.  Assume  as  boundary  conditions  a  distribution  of  blowing 
velocity  q  (x,0)  along  the  plate  and  a  tangential  velocity  q  (x,0)  =  0. 

y  ^ 

The  latter  condition  is  not  necessary  but  is  believed  to  correspond 
to  the  correct  limit  of  the  viscous  flow  equations  unless  special  pre¬ 
cautions  are  taken  with  the  blowing.  These  boundary  conditions  intro¬ 
duce  rotation  into  the  injection  layer.  The  additional  boundary  condi¬ 
tion  can  be  accommodated,  since  the  location  of  the  dividing  streamline 
is  not  known  in  advance  but  must  be  found  as  part  of  the  solution  to 
the  problem.  The  conditions  to  be  satisfied  across  the  separating 
streamline  are  continuity  of  pressure  and  direction  of  flow. 

Let  the  dividing  streamline  be  represented  by 

y  =  6S(x)  (2.1) 

where  S(0)  =  0,  S(l)  =  1,  and  6,  the  basic  small  parameter  of  the  prob¬ 
lem,  is  the  thickness  of  the  blown  layer  at  a  unit  distance  downstream 
of  the  nose.  The  thickness  of  the  layer  tends  to  zero  as  6  -*  0,  cor¬ 
responding  to  the  disappearance  of  both  blowing  and  the  pressure  and 
tangential  velocity  perturbations.  In  order  to  obtain  approximate 
equations  valid  within  the  thin  layer,  distances  normal  to  the  wall 
must  be  measured  in  terms  of  the  thickness  of  the  layer.  That  is,  a 
coordinate 


y 


(2.2) 


is  held  fixed  as  6  -*  0 .  Assume  now  that  the  flow  in  the  thin  layer  is 
incompressible  (more  precise  conditions  for  the  validity  of  this  assump¬ 
tion  will  be  given  later) .  Thus  the  following  forms  of  asymptotic  ex¬ 
pansions  in  6  can  be  assumed  for  the  flow  quantities  in  the  thin  layer: 
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(x,y;fi)  =  3(S)v(x,y)  + 


(2.4) 


P(x,y;6)  -  Pro  „ 

- 2 -  =  e<6)p(x>y)  + 

P  U 


(2.5) 


The  continuity  equation  is  nontrivial  for  (3  =  6a,  and  the  tangential 

2 

momentum  is  balanced  between  inertia  and  pressure  gradient  for  e  =  a  ; 


3  = 

e  = 


(2.6) 


The  resulting  approximate  equations  are 


continuity: 


by 

dy 


0 


(2.7) 


du  ,  du  dp 

x-momentum:  u  -r-  +  v  —  =  -a  -5“ 

dx  dx 

dy 


(2.8) 


y- momentum:  0 


(2.9) 


where 

a  =  density  ratio  (of  order  one)  =  Pc0/Pw. 

Thus  we  obtain  essentially  the  inviscid  version  of  the  Prandtl 
boundary- layer  equations.  Hence  the  results  above  still  apply  to  a 
curved  surface  with  an  arc- length  x  and  a  normal  distance  y,  so  long 
as  the  radius  of  curvature  is  larger  than  0(6)  in  characteristic  units. 

A  consequence  is  that  the  pressure  in  the  layer  p  is  only  a  func¬ 


tion  of  x: 
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P  =  P(x) 


(2.10) 


and  we  obtain  the  further  simplified  equations 


du 

dx 


+ 


Sv 

3y 


0 


(2.11) 


u 


3u  .  3u 

"S — H  v  — 

ox 

dy 


(2.12) 


We  can  now  estimate  the  Mach  number  in  the  inner  layer  by  using 


uV(6)pra 

YCTpJl  +  ..  .] 


mV  (§) 


(2.13) 


We  see  that  if  the  external  Mach  number  is  fixed  as  6  -»  0,  then  the 
flow  is  incompressible,  although  compressibility  must  be  taken  into 
account  in  the  hypersonic  case  (M  -*“>).  We  can  also  estimate  the 
neglected  viscous-stress  terms  in  Eq.  (2.8)  by  calculating  the  ratio 


-.2 

o  q 

jo_ _ _2S 

P  a  2 
w  oy 

3q 

x 

^x  dx 


V 


<x(6)6  l 


(6  a) (Re) 


(2.14) 


We  assume  that  the  Reynolds  number  based  on  the  characteristic  length 

(■£  =  1)  is  large  enough  to  make  the  ratio  in  Eq.  (2.14)  much  less  than 

2 

one.  For  the  theory  to  apply,  Re  must  be  regarded  as  Re(5)  »  1/(6  a). 
Thus  for  a  given  Reynolds  number  the  theory  is  consistent  when 


—  «  ad2  «  l 
Re 


(2.15) 


That  is,  the  layer  must  be  thin,  but  not  too  thin. 
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The  system  of  Eqs .  (2.11)  and  (2.12)  is  now  to  be  solved  with  the 
boundary  conditions 


? 


v(x,0)  =  vw(x),  u(x,0)  =  0 


(2.16) 


It  is  convenient  at  first  to  regard  P(x)  as  known  and  to  see  if  some 
relation  can  be  found  among  v^(x) ,  P(x),  and  S (x) .  This  is  most  easily 
done  if  we  introduce  a  new  coordinate  that  is  constant  along  a  stream¬ 
line  in  the  blown  layer.  Let  the  stream  function  i|i(x,y)  be  defined  by 


u(x,y)  =  ,  v(x,y)  =  - 

dy 


(2.17) 


Then  on  the  surface  (y  =  0) 


v  (x) 
wv  ' 


(xj0) 


or 


^(x)  =  Kx,0)  = 


-I  vK«)d{ 
0 


(2.18) 


* . 


Let  x  (\|i)  be  the  x-coordinate  where  a  streamline  enters  the  flow  at 

y  =  0;  that  is,  the  inverse  function  of  Eq.  (2.18)  (see  Fig.  2). 

Using  v  to  represent  the  velocity  at  the  plate  we  have 
w 


* 


dx 


dty 


or 


}  dY 

x  <*>  =  -J  TTf) 
o  w 


> 


(2.19) 
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The  problem  can  be  formulated  in  the  (x  ,x)  plane  (Fig.  2)  where  the 

„  * 
surface  y  =  0  is  x  =  x,  the  separating  streamline  is  x  =0,  and 

* 

0  ^  x  ^  x. 

We  apply  the  following  transformations  from  (x,y)  to  (x,x  ): 


d  ^  d  dx*  5 

c)x  *5Z  dx 


a  dx*  a 


v(x,x  )  d 

v  (x*)  dx* 
w'  ' 


(2.20) 


d  dx*  d 

dy  dy  dx 


u(x,x  )  d 
*  .  * 
v  (x  )  dx 
wv 


(2.21) 


u 


d 


d 


+  v  — 
dy 


*  d 
u(x,x  )  ^ 


(2.22) 


Thus  the  tangential  momentum  equation  (2.12)  becomes,  in  the  new  coordi¬ 
nate, 


du  dP 

U  dx  =  "a  dx 


(2.23) 


This  can  be  integrated  to  give  the  approximate  version  of  the  Bernoulli 
equation  if  the  no-slip  condition 

r"  ) 

u(x*,x*)  =  0  (2.24) 

is  taken  into  account.  This  approximate  Bernoulli  equation  is: 

I  u2(x,x*)  +  CTP(x)  =  CTP(x*)  (2.25) 

The  horizontal  component  of  the  flow  is 


u(x,x*)  =  (2a)1/2(p(x*)  -  p(x))1/2 


(2.26) 
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This  result  implies  that  the  solution  is  valid  only  for  a  favorable 

.  * 

pressure  gradient  with  P(x  )  >  P(x). 

Next  the  shape  of  the  dividing  streamline  is  calculated  by  inte¬ 
gration  across  the  streamlines. 

For  x  =  constant, 


dy 


d\|r  _  di|f  dx 

u  ,  *  u 

dx 


dx 


* 

u(x,x  ) 


(2.27) 


Thus,  integrating  Eq.  (2.27),  with  the  boundary  condition  that  the 
dividing  streamline  originates  at  the  nose  [y  =  S(x)  for  x  =  0] ,  we 
obtain 


y(x,x  )  =  s(x)  -  J 


X*  V  (!) 


w ' 


u(x,|) 


d§ 


(2.28) 


Using  the  Bernoulli  equation  (2.26),  we  replace  the  tangential  velocity 
in  the  last  expression  by  its  equivalent  in  terms  of  pressure; 


* 


y(x,x  )  =  s(x) 


X 

S 


V  (§) 

wv  ' 


(2ct)1/20  (p(|)  -  P(x) 


1/2 


d?> 


(2.29) 


Finally,  this  last  expression  can  be  evaluated  on  the  surface  where 
*  ~ 

x  =  x  and  y  =  0  to  provide  a  relationship  among  S(x),  v^(x) ,  and  P(x): 


S(x) 


J 


0 


Vs> 

(p(|)  -  P(x))1/2 


d| 


(2.30) 


In  general,  the  relation  between  the  pressure  and  the  shape  of  the 
dividing  streamline  depends  on  the  nature  of  the  external  flow.  It 
is  therefore  convenient  to  express  vw(x)  as  an  integral  over  P,S,  since 
Eq.  (2.30)  is  an  integral  equation  of  Abel  type  for  vw(5)  (see  Appendix 
for  details) : 
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v 

w 


(X) 


-  1  (2.)1/2  £ 

TT  v  '  dx 


llSSL 


o  (p(S)  -  poo) 


1/2 


d§ 


(2.31) 


Equation  (2.31)  is  the  basic  answer  to  our  problem,  and  provides  the 

desired  information  in  an  inverse  form.  The  blowing  distribution  v  (x) 

w 

can  be  found  for  a  given  P(x),  S(x). 

We  shall  now  apply  this  theory  concretely  by  specifying  the  outer 
flow  and  its  perturbations.  First,  we  can  fix  the  order  of  magnitude 
a,|3,e  of  the  perturbation  by  constructing  a  model  of  the  outer  flow. 

For  example,  consider  flow  past  a  flat  plate  and  conditions  for  which 
linearized  subsonic  or  supersonic  theory  is  valid;  that  is, 


«  1 


(2.32) 


Then  it  is  well  known  that  the  pressure  perturbations  in  the  outer  flow 
are  of  the  same  order  as  the  body  thickness  or  shape;  that  is, 

e(6)  =  6 

1  Qf(6)  =  61/2  ,  (2.33) 

p(6)  «  63/2 

The  last  equation  relates  d  to  the  order  of  magnitude  of  the  blowing 
velocity: 


6 


(^) 


(2.34) 


To  ensure  the  self-consistency  of  the  theory,  the  following  more 
precise  estimates  of  Reynolds  number  and  blowing  velocity  can  now  be 
obtained  from  Eq.  (2.15): 
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or 


«  6  «  1 


« 


qy(x,o) 

u 


«  i 


If  this  theory  is  to  supersede  classical  boundary-layer  theory 

1/2 

(for  a  =  1),  it  is  necessary  that  q  /U  >  1/ (Re)  .  This  condition 

ensures  that  the  thickness  of  the  blown  layer,  6,  is  greater  than  a 

1/2 

characteristic  thickness  of  the  shear  layer,  l/(Re) 

Thus,  for  linearized  external  flow,  the  basic  expansion  [(Eqs. 
(2.3)  through  (2.5)]  becomes 


4x  1/2  ~ 

-y-  (x,y;6)  =  6  u(x,y)  + 


^  (x,y;6)  =  63/2  v(x,£)  + 


p(x,y;S)  -  pro 

- ^ -  =  6p(x)  +  ... 


P  XT 
r  00 


(2.35) 

(2.36) 

(2.37) 


The  simplest  results  are  obtained  for  linearized  supersonic  (Ma  ,  >  1) 
external  flow,  where  the  following  relation  between  pressure  perturbation 
and  local  slope  is  valid: 


P(x) 


(2.38) 


Thus,  the  blowing  velocity  and  shape  of  the  dividing  streamline 
can  be  expressed  in  terms  of  a  given  pressure  distribution  from  Eq. 
(2.31)  : 
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v 


w 


«  -  -  *  (o4  -  I I 


PHI 


0  (*(§)  -  P(x)) 


Y72  (2.39) 


/  o  \l/2  X 

S(x)  =  (Mro  -  1 J  J*  P(S)d§ 

0 


(2.40) 


We  next  apply  the  foregoing  method  to  some  specific  examples. 

Example  1  --  Ogive  Nose 

Let  the  pressure  distribution  be  given  by 


P(x)  =  Sq  -  a^x  0  <  x  <  1 


(2.41) 


Then  Eq.  (2.39)  gives 


1/2  x  a  -  a  § 


i  /  9  \  '  n  i  = 

4  -  » )  I  /  °  -.V\i72 

0  Vax(x  -  l) J 


d5 


,1/2 


-  -i  (2»od  - 1))  <V1/2  K*1/2  - !  V3/2}-  <2-«> 


and  Eq.  (2.40)  gives 


s(x)  =  (*£  -  0  (v  "  \  aix2} 


(2.43) 


The  normalization  S(l)  =  1  gives  the  following  relation  between 
a^  and  a^; 


(2.44) 
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However,  the  condition  of  favorable  pressure  gradient  requires  that 

a  >  0 


so  that  the  solution  is  only  valid  for 


(2.45) 


It  is  seen  from  this  example  that  an  ogive- like  layer  with  a  con- 

2  1/2 

stant  slope  6 (M  -  1)  a  at  the  nose  is  produced  by  blowing  which 
1/2°°  0 

causes  ~  (x)  to  vanish  as  x  -*  0 ;  however,  it  is  the  curvature  of 
the  ogival  layer  near  the  nose  which  determines  the  required  injection 

•j— ^  j  _  ~j 

distribution  as  x  -*  0.  For  example,  for  S  to  be  Ax  -  Bx  ,  and  P  to 

/  9 

be  A(n  +  l)Bx  ,  it  is  necessary  that  vw  ~  x  as  x  -*  0.  The  injection 
velocity  must  vanish  as  x  -*  0  in  order  to  maintain  the  finite  slope  of 
the  injection  layer. 


Example  2  --  Power-Law  Nose 
Let 


„ ,  N  1  -  n  -n 

P(x)  =  - T-7T7  x 


(m2  -  1 

\  CO 


1/2 


(2.46) 


S(x)  =>  x1  n  n  <  1 


(2.47) 


Then 


x 

J 


r-n  ,  .1/2  r(i  -  1)  1 

i  dc  (tt)  \n  2/  1- 


n/2 


0  (5-n  . 


1/2 


n 


0  <  n  <  2  (2.48) 
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so  that 


(2.49) 


There  is  a  finite  rate  of  mass  addition  only  for  n  <  2/3.  The  often- 

-1/2  -1/3  2/3 

used  case  where  vw  ~  x  corresponds  here  to  p  ~  x  ,  S  ~  x 

Note  however,  that  because  P  -*  <»  as  x  -»  0,  the  Bernoulli  integral  (2.26) 
indicates  that  u  -*  ro  on  the  dividing  streamline.  The  solution  is  prob¬ 
ably  not  valid  near  the  dividing  streamline  for  these  cases,  and  the 
calculation  of  a  free  shear  layer  would  require  further  discussion. 

The  order  of  magnitude  of  the  jump  in  tangential  velocity  at  the 

separation  streamline  can  be  calculated.  Since  the  external  velocity 

1/2 

is  U  +  0(6)  and  the  tangential  velocity  in  the  blown  layer  is  0(6)  , 

_  * 

the  jump  in  tangential  velocity  LqJ  is  expressed  in  terms  of  u(x,x  )  as 

-5  [q]Tang  =  1  -  (6)1/2u(x,0)  =  1  -  (26ct)1/2[P(0)  -  P(x)]  (2.50) 

Within  the  context  of  this  thin- layer  theory,  the  vorticity  is  always 
finite  and  is  introduced  at  the  wall  by  the  stagnation  pressure  gradi¬ 
ent  of  the  injectant;  the  skin  friction  coefficient  is 


1/2 


7T72I 

6  Re. 


w 


dp 

dx 
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III.  OTHER  EXAMPLES  OF  INCOMPRESSIBLE  INJECTION  LAYERS 


In  this  section  we  consider  two  further  examples  of  inviscid  flows 
with  mass  transfer  at  the  surface,  which  are  characterized,  like  the 
flat  plate  in  linearized  supersonic  flow,  by  incompressible  flow  in 
the  thin  injection  layer. 

HYPERSONIC  WEDGE  FLOW  WITH  DISTRIBUTED  BLOWING 

The  simple  theory  presented  in  Section  II  can  be  readily  extended 
to  the  problem  of  supersonic  flow  past  an  inclined  flat  plate  (wedge), 
with  distributed  mass  injection  along  the  surface  of  the  wedge.  We  con¬ 
sider  only  the  case  where  the  injection-layer  thickness  6  is  small  com¬ 
pared  to  the  original  wedge  height  9. 

According  to  the  theory  of  supersonic  flow  past  a  slightly  per¬ 
turbed  wedge,  flow  quantities  in  the  outer  flow  can  be  expressed  as 


qx(x,y) 

-  U2[l  +  eu(x,y)] 

(3.1) 

qy(x,y) 

=  U2ev(x,y) 

(3.2) 

p(y) 

=  P2  +  ep2U2p(x,y)  +  . . . 

(3.3) 

where  e  is  a  small  parameter  determined  by  the  magnitude  of  the  pertur¬ 
bation,  and  the  subscript  "2"  refers  to  conditions  behind  the  original 
shock  wave. 

The  condition  of  tangential  flow  at  the  outer  edge  of  the  injectant 
layer  is  replaced  by  a  boundary  condition  on  the  original  wedge  surface, 
and  the  perturbed  shock  jump  conditions  are  given  along  the  original  un¬ 
perturbed  shock  wave  of  the  wedge.  Reference  10  contains  a  complete 
discussion  of  this  theory. 

A  solution  of  the  wave  equation  with  appropriate  boundary  condi¬ 
tions  results  in  the  following  equation  for  the  pressure; 


P 


+ 


o1/2 


+22  Xny 
n=l 


(3.4) 
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where 

X  =  reflection  coefficient 

k  =  ratio  of  the  x-coordinate  where  a  wave  originates  near  the 

wedge  to  the  x-coordinate  where  the  wave  intersects  the  wedge 
y(x)  =  dimensional  perturbation  measured  from  the  wedge  surface. 


If  this  result  is  specialized  to  the  hypersonic  limit,  Mro  -*  », 

9  ,  0  M  9  ,  -*  =°,  then 

wedge  00  wedge 


X  -*  A.  = 

00 


.  _ ) 

\2(Y  -  1)J 


1/2 


1  + 


N2  (Y  -  1) 


1/2 


(3.5a) 


k  -*  k 


a/2 


1  + 


fy  -  i 

\  2Y 


1/2 


(3.5b) 


Equation  (3.4)  becomes 


p(x,6s(x) J  =  (Y  +  1)9 


p  u2 

1  00  CO 


1  + 


6/  2Y 

9  W  -  1 


1/2 


S  /(x)  +  2  L  XnS''(knx) 
n=l 


(3.6) 


where  the  shape  of  the  bounding  streamline  is  Y(x)  =  0(6/9)S(x),  and 
6/9  «  1. 


Inner  Flow 

The  injected  fluid  is  assumed  to  leave  the  surface  with  a  constant 
temperature  T  .  This  assumption  is  not  essential  to  the  analysis,  how¬ 
ever,  and  the  theory  remains  valid  for  moderate  variations  of  wall  tem¬ 
perature. 

We  introduce  a  new  coordinate  y  =  y/6.  Consistent  orders  of  mag¬ 
nitude  for  flow  variables  in  the  injection  layer  are  as  follows: 
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qx(x,y) 


u(x,y)  +  .  .  . 


qy(x>y) 


^w 


1/2 


3/2 


(3.7a) 


(3.7b) 


Y  +  1  ^ca^co  2~  ~ 

p(x,y)  =  y  -  -p  'nr-  6  P(x,y)  +  ...  (3.7c) 


P(x,y)  =  ^P^e2  {l  +  |Pg(x,y)+  ...}  (3 •  7d) 


where  is  the  universal  gas  constant  and  u  is  the  molecular  weight 
0  w 

of  the  injectant.  The  equations  of  motion  then  become 


x- momentum: 


+  v 


= 

dy 


(3.8) 


continuity: 


5(uP)  +  Q(vP) 
dx 

dy 


(3.9) 


dps 

y- momentum:  -  =  0 

dy 


(3.10) 


If  we  assume  a  perfect  gas  with  an  isentropic  exponent  Y  ,  then  the 
entropy  equation  can  be  written 


4;  =  £(« 

P 


If  the  wall  temperature  is  constant  and  the  injection  Mach  number  is 
small,  then  the  flow  is  incompressible,  since  p  -  1  +  0(6/9).  The 
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boundary  conditions  at  the  wall  become  u(x,0)  =  0,  v(x,0)  =  v^(x) . 
The  boundary  conditions  on  the  dividing  streamline,  y  =  S(x),  are 


v(x,S(x)) 

u(x,S(x>) 


S'(x) 


(3.11) 


and 


P 


S 


S  '(x)  +  2  E  XnS  '(knx) 
n=l 


(3.12) 


As  in  Section  II,  the  equations  of  motion  and  entropy  reduce  to  those 
of  an  inviscid,  incompressible  boundary  layer.  The  density  is  constant, 
and  there  is  no  pressure  gradient  across  the  layer. 

In  terms  of  the  mass  flow  ratio  per  unit  area,  m  =  (pq  )  /p  U  , 

*■  3  \  i  1  CO  CO 

the  requirement  that  6/8  «  1  results  in  the  following  estimate: 


or 


6_ 

9 


m 


££)  ]  «> 


1/2  2/3 


0' 

00 


Ui 
oo 


or 
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For  a  fixed  Mach  number  M^,  a  fixed  mass  flow  m,  a  fixed  wedge  angle  9, 
and  a  fixed  wall- to- free- stream  temperature  ratio  T^/Tm,  this  implies 
that 


6 

0 


1/2  1/3 

] 


Thus,  the  injection  of  a  light  gas  into  a  heavy  one  produces  greater 
induced  pressure  effects  than  does  the  injection  of  a  heavy  gas  into 
a  light  one.  For  example,  pressure  perturbations  due  to  the  injection 
of  helium  into  air  are  predicted  to  be  about  twice  those  for  air  into 
air . 

The  flow  in  the  injection  layer  is  governed  by  the  same  equations 
as  those  given  earlier,  and  we  obtain  the  simple  Bernoulli  equation: 


y-  (x,x*)  +  pg(x,x*)  =  pg(x*)  (3.13) 

where  x^  is  defined  in  Eq.  (2.19).  The  height  of  the  dividing  stream¬ 
line  is  again 


S(x) 


,1/2 


x 

J\ 

0  "  ?s 


(x)J 


1/2 


d§ 


(3.14) 


where 


1/2  CO 

ps(x)  =  (y  i)  {s'(x)  +  2  E^^S^^x)}  (3.15) 

As  before,  an  inverse  representation  of  the  solution  is  possible:  S(x) 
is  assumed,  the  resulting  pressure  distribution  is  calculated  from  Eq. 
(3.15),  and  Eq.  (3.14)  then  becomes  an  Abel  integral  equation  for  vw(x) . 
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Example 


If 


_  /  y  .  ^-1  n 

S  (x)  =  Sax 
n  n 
n=0 


(3.16) 


then 


Ps(x,0) 


(  2  )1/2 

N' 

V 

1  + 

A  knl 
00  00 

n 

o  v 

\Y  -  1  / 

n=0 

.1  - 

A  kn 

00  00  J 

a  a. 

n 

(3.17) 


If  the  blowing  results  in  an  ogive- like  injection  zone. 


n  x 

S  —  a^x  -  a 2  2 


S  =  a1  -  a^x 


/  where  a^,  > 


(3.18) 


The  pressure  formula  gives 


f  2  )1/2 

l1  + 

A 

CO 

1  +  A  k 

00  CO 

\Y  -  1/ 

1 1  - 

T~  ai  " 

1  -  A  k  a2 

,  vw(x)  is 

then 

CO 

00  00 

(3.19) 


v  =  + 

w 


tt  LVY 


„  \  1/2  1  +  A  k 

2  \  00  CO 


1/2 


=—t  X  - - —  X  an  X  2 

-  1/  1  -  A  k  2 


x  (a  -  a  A) 

J  . X7_  | .  (3.20) 


or 


v  (x) 
wN  ' 


TT  L  VY  -  1 


2  ^2  1  +  KK 

1  -  X  k  x  a2  x  2 
00  00 


^2[x1/2 


X  a  x  2  -  X3/2  X  a2  X 
(3.21) 
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Since  S(l)  =  1,  then 


a 


1 


=  1 


(3.22) 


or 


a2  =  2{ax  -  l}  (3.23) 

HYPERSONIC  NEWTONIAN  FLOW  PAST  A  POROUS  CONE 

Another  example  of  the  "blowhard"  problem  is  provided  by  super¬ 
sonic  or  hypersonic  flow  past  a  porous  cone  when  the  injection  rates 
are  high.  We  restrict  our  attention  to  the  study  of  hypersonic  flow 
past  a  porous  slender  cone  of  angle  0,  and  use  the  Newtonian  formula 
for  the  calculation  of  the  pressure  distribution  on  a  perturbed  cone. 

As  before,  the  "blown"  layer  is  assumed  to  be  thin  relative  to  the 
thickness  of  the  original  cone. 

Newtonian  slender- body  flow  corresponds  to  the  limits  0  -*  0, 

2  2 

Y  -*  1,  and  1/M  0  (Y  -  1)  -»  0.  The  Newtonian  formula,  including  impact 
and  centrifugal  force  terms,  is 


P 


(3.24) 


where  the  profile  shape  is  r  =  r(x)  and  primes  indicate  d/dx. 

If  the  injectant  layer  height  at  x  =  1  is  6,  then 

r  =  9[x  +  (§)  S(x)]  (3.25) 


where  6/0  «  1.  The  initial  height  and  the  normalization  condition 
are  given  by: 
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S(0)  =  0,  S(l)  =  1 


(3.26) 


This  pressure  formula  then  becomes 


P 


Y1  ®2  {l  +  2  \  S  '(x)  +  (!)  x  (x)  +  0 


(3.27) 


The  following  orders  of  magnitude  are  appropriate  for  the  flow 
variables  in  the  inner  layer: 


q 


x 


q 


y 


v(x,y) 


(3.28a) 


(3.28b) 


P 


2  2 
p  U  0 

1  OO  00 


[’ 


+  9  Ps(x,y) 


+ 


(3.28c) 


p  U  92 

■  CO  CO  _  ^ 

P  =  ~2r"t -  P(x>y) 

0  w 


(3 . 28d) 


The  equations  of  motion  become,  in  boundary -layer  coordinates. 


_d_ 

dx 


du  ,  du 

u u  +  ■ 

[r0(x>u]  +  z:  [ro(x)v]  ■ 

dy 

P  = 


dPo 

(3.29) 

dxp 

0 

(3.30) 

1 

(3.31) 
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dp£ 

Sy 


=  0 


(3.32) 


where  y  =  y/6  and  r^(x)  =  x9. 

By  introducing  a  stream  function  defined  by 


-  %<*>•  *x 


-vr0(x) 


(3.33) 


we  obtain  the  simplified  Bernoulli  equation  of  Section  II: 


4  (t  (xJ)  +  =  0  (3-34> 

We  also  introduce  the  variable  x  ,  which  marks  the  location  of  the 
origin  of  the  streamline  t  on  the  surface  of  the  cone: 


,  *  * 
rQ(x  )dx  = 


d{' 


(3.35) 


The  height  above  the  cone  surface,  y(Xj'lf),  can  be  obtained  from 


or 


ro<x)dy  '  TT 


(3.36) 


r0(x)y(x,x  ) 


*  r0(l)vw(l) 

rQ(x)S(x)  -  f  - rrr  .a-  d? 


0 


u(§,x) 


(3.37) 


Introducing  the  Bernoulli  equation,  Eq.  (3.34),  we  obtain  the  height 
S(x): 
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r0(x)S(x)  = 


x  rn(§)v(|) 


(2) 


1  }  0 '  '  w  ^  y 

1/2  o  (psd)  -  ps(x>)1/2 


d§ 


(3.38) 


We  now  consider  a  solution  to  the  inverse  problem  of  finding  the 
blowing  distribution  which  leads  to  a  prescribed  S(x). 

From  the  Appendix, 


r0(x)Vw(x) 


(2) 


1/2 


i  dP  }  dfVro(§)s(5); 

■  "  dl  l  (psCS>  -  PSWW/2 


d§ 


(3.39) 


An  especially  simple  example  of  the  application  of  the  theory 
corresponds  to 


P<x)  =  a0  "  aLx 


(3.40) 


then 


.  30  31  2 
S(x)  =  —  x  -  —  x 


V  (x) 
WN  ' 


,1/2 


<*!  X  2)*'"  i  {a0x‘'-  X  5  -  a.x"-  X 


.1/2  „  4 
3 


3/2 


-I 

25J 


(3.41) 


5 

2  al 


5  >  0 


Since  S(l)  =  1,  then 
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IV.  CCMPRESSIBLE  INNER  LAYERS 


If  the  limit  of  the  gas  density  as  6  ->  0  is  not  a  constant,  we 
must  consider  the  effect  of  compressible  inner  layers.  Examples  of 
these  are  surfaces  which,  in  the  absence  of  injection,  have  favorable 
pressure  gradients,  or  situations  in  which  the  hypersonic  parameter 
MS  -*  ».  In  this  section  we  derive  the  simplified  inviscid  equations 
of  motion  corresponding  to  a  compressible  inner  layer.  The  free  stream 
is  assumed  to  be  hypersonic. 

The  following  consistent  set  of  expansions  is  assumed  valid  as 

~  2  2 
6  -*  0,  with  y  fixed  (it  is  assumed  that  Mra  -*  00 ,  11^6(6)  -*  »,  and  6  -*  0)  : 


qx  =  (nr2)  v-v 

w 

R  T  1/2 

qy  =  (nr2)  6v<x’y>  (4-2> 

w 

T 

Y-  =  t  (x,y)  (4.3) 

w 

P(x,y)  =  P0OU2e(6)p(x,y)  (4.4) 


P 


pjsI 

RnT 
0  w 


w 


e(6)P(x,y) 


(4.5) 


Y~  =  t(x,y)  (4.6) 

w 


where  x  and  y  are  boundary -layer  coordinates  along  and  normal  to  the 
surface.  With  these  expansions,  the  momentum  equation  balances  inertial 
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and  pressure-gradient  forces,  the  continuity  equation  is  preserved, 
and  the  normal  pressure  gradient  dp/dy(x,y)  0  as  6  approaches  zero. 
The  momentum  equation  becomes 


du  ,  du  1  dp 

U  3 - h  V  3—  =  -  —  33 

ox  oy  _  dx 

J  P 


(4.7) 


The  continuity  equation  remains 


3  d  _ 

^  (pu)  +  —  (pv) 

dy 


=  0 


(4.8) 


and  the  energy  equation  becomes 


ato  dto 

U  3"-  ■  +  V  3 -  =  0 

dx  dy 


(4.9) 


The  total  enthalpy,  in  nondimensional  form,  is  then 


tQ(x,y)  =  t(x,y)  +  u2(x,y) 


(4.10) 


and  the  equation  of  state  becomes 


p  = 

t(x,y) 


(4.11) 


Von  Mises  variables  are  introduced,  replacing  x,y  by  x,t.  The  trans¬ 
formation  formulas  are 


'l'-  =  pu 

y 


(4.12a) 


♦ 


x 


-pv 
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The  momentum  equation  becomes 


u  +  =  0 
dX  pSx 


(4.12b) 


(4.13) 


and  the  energy  equation  becomes 


u  =  0  (4.14) 

The  integral  of  the  total  energy  equation  is  then  t^x,'!')  =  t^(t) ,  where 
t0(t)  is  a  known  function  of  t  corresponding  to  the  wall  temperature. 

The  momentum  equation  can  be  integrated,  giving: 


V*)  ‘  T  (x’^  =  K(^  x  p<x»+> 


Y-l 

Y 


(4.15) 


Changing  variables  to  x,x  ,  where  x  marks  the  location  on  the 
surface  where  the  t  streamline  enters  the  flow,  we  obtain: 


ft  ^ 

p(x  ,x  )vw(x  )dx  "  =  -dt 


(4.16) 


*  *  * 

At  a  point  on  the  surface  where  x  =  x  ,  u  =  0  and  p  =  p(x  ,x  ) . 


Thus 


Y-l 


u2(x,x*)  =  ^0(x*)X*)  X  {l  -  [  } 


Ll 

ft 

‘P(x",x  )* 


(4.17) 
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In  order  to  locate  the  point  x,y  corresponding  to  x,x  ,  we  use 


*  *  „  *  *  , 

Ip.Ql dx ' 

"k  k 

p(x  ,x)u(x  ,x) 


(4.18) 


Noting  that 


1/Y,  *  •, 
P  <2 L  »x) 

-/  * 
p(x  ,x) 


1/Y,  *  * 
P  (x  »x  ) 

p(x  ,x  ) 


(4.19) 


*  *  1/Y  .  * 

i~  /  *N  p(x  ,x  )  dx 

dy  =  -v  (x  )  rA-  .  ? — ^  - t - 

w  /  *  \  \ 
P (x  ,x)  u(x  ,x) 


(4.20) 


Equation  (4.20)  can  be  integrated  to  yield 


^-il/Y 


y(x*,x)  =  S(x)  -  J  vw(§)[-^|^-_ 


vu(§»x) 


(4.21) 


Since  y(x,x)  =  0, 


■»  ■ 

\Y  -  1/ 


«b«> 1  - 


(4.22) 
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Example  --  Power-Law  Nose,  e(6)  =  6 
Let 


S  =  xn,  n  <  1  (4.23) 

For  a  power- law  effective  nose  shape,  hypersonic  small-disturbance  theory 
yields 


P 


A  X  x2^1) 
n 


(4.24) 


The  injection  distribution  is  v =  v^x111;  Eq.  (4.22)  is  used  to  determine 

v  and  m  as  follows: 
m 


H-)U\  x" 

.7  -  1/  0* 


=  v 


m 


x  r 


II  Y 

x_ 


d§ 


-  (?) 


2 (n-1)  X 


Y 


-TT72  <4'25> 


or 


1/2  1 
2Y  V  .  v  n  v  nri-1  f* 

- -  )  tn  X  x  =  v  X  x 

-1/0  m  J 

0 


mf  -  (n-1) 


i  -  (X)-2*-1) x  ¥ 


,1/2 


dA 


(4.26) 


where  X  =  §/x.  Therefore 


n  =  1  +  m 


v 

m 


(4.27) 


(4.28) 
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where 

A  =  m  +  |  (n  -  1) 

B  =  -2(n  -  1)  X 

We  note  that 

(4.29a) 


(4.29b) 

where  T  is  the  gamma  function. 


XAdX 


j. 

s 

0  (l  -  X 


B\ 


1/2 


=  ir 

B  J 


dx  X  x 


Afl-B 
B 


1  -  x 


1/2 


r/A  +  l\  r(l\ 

1  V  B  /  \2  J 


B  n/A  +  1  1 

1  B  +  2 


Example  --  Two-Dimensional  Stagnation  Point 

Although  this  is  not  strictly  a  blowhard  problem,  it  can  be  treated 
by  the  thin- layer  inviscid  theory  as  follows: 


s(6)  -  1 


P 


PQ[1  -  ^2]  +  ... 


(4.30) 


v  =  constant 
w 


(4.31) 


S(x) 
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S(x) 


v 

w 


(2) l/2(tQa) 1/2 


w 


(2)1/2(t0a)1/2 


(4.32) 


(4.33) 


(4.34) 
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V.  COMPRESSIBLE  AXISYMMETRIC  FLOW 


The  rotationally  symmetric  version  of  the  compressible  blowhard 
equations  for  a  hypersonic  free  stream  is  derived  in  this  section.  A 
boundary- layer  coordinate  system  is  used  and  an  additional  parameter 
R/6  is  introduced  (see  Fig.  3). 


r  =  y  cos  0  +  Rr^(x) 


^  ^  R 

r  =  y  cos  0  +  j  rQ(x) 


(5.1) 


(5.2) 


w 


tQ(x,y) 


(5.3) 


P(x,y)  =  ProU^e(6,|)p(x5y) 


(5.4) 


f~  =  t(x,y)  (5.5) 

w 


The  inviscid  balance  between  inertia  and  pressure  gradient  becomes 


du  ,  du 
u  -r — 1  v  — 
ox  -n~ 
dy 


1  d£ 

p  dx 


(5.6) 


The  continuity  equation  becomes 
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Fig.  3  --  Coordinate  system  for  axisymmetric  flow. 


-37- 


(urp)  +  (vrp)~  =  0 

x  y 


(5.7) 


and  the  energy  equation  becomes 


ut„  +  vt.~  =  0 

Ox  Oy 


(5.8) 


The  equation  of  state  is  p  =  p/t.  The  boundary  conditions  are; 


y  =  0,  v  =  vw(x),  u  =  0,  tn  =  tnii  (5.9) 


0  Ow 


y  =  S, 


v  dS 
u  ~  dx 


Introducing  the  stream-function  variables 


pur  =  l|l; 


(5.10) 


pvr  =  - t 


x 


we  obtain  the  compressible  Bernoulli  equation,  which  gives  us 

Y-l 


2  *  *  y 

u  (x  ,  x)  =  2tQ(x  )  X  — — — 


Y 


i  _  p.fcjX-L- 
*  *. 
p(x  ,x  ) 


(5.11) 


where 


* 


dx 


d'l' 


_  *  *  ,  * 
P(x  ,x  )  X  vwrQ(x  ) 


(5.12) 


In  order  to  locate  the  point  (x  ,x) ,  we  must  obtain  the  y  coordinate 
from 
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rdy  = 


dt 

Pu 


(5.13) 


or 


2“  cos  8(x)  +  |  rQ(x)y 


cos  9(x)  +  r^S  X  S 


r0vwP(§V) 

“  J  — -  d? 

p(§,x)u(§,x) 


(5.14) 


The  edge  of  the  "blown  layer"  is  determined  from 


§1  cos  9  +  £  r  S  =  *  }  r0Vwdg  (2 ILii) 
2  S  0  6  J  u(§,x)  Vp(§,x)y 


,1/Y 


(5.15) 


For  practical  applications,  it  is  anticipated  that  R/6  ~0(1).  Further 
analysis,  involving  the  introduction  of  a  pressure  formula  for  the 
outer  flow,  leads  to  the  formulation  of  an  inverse  problem  like  those 
in  the  preceding  sections. 
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VI.  CONCLUSIONS 


A  simple  theory  has  been  proposed  for  the  analysis  of  hypersonic 
or  supersonic  flows  over  cones,  flat  plates,  and  wedges,  where  there 
is  continuous  mass  injection  from  the  surface.  The  mass  transfer  rates 
are  assumed  to  be  low  enough  to  preserve  the  slenderness  of  the  effec¬ 
tive  body  produced  by  the  injection  but  high  enough  to  lift  the  boundary 
layer  off  the  wall.  Since  a  favorable  pressure  gradient  is  necessary 
to  drive  the  flow  in  the  injection  layer,  the  theory  applies  only  to 
those  injection  distributions  that  cause  the  pressure  distribution  to 
fall.  Hass  transfer  must  be  provided  externally;  among  the  possible 
methods  are:  surface  ablation  caused  by  externally  applied  thermal 
radiation,  transpiration  through  porous  walls,  or  a  reaction  occurring 
at  the  boundary,  Hass  transfer  caused  by  convective  heat  transfer  is 
excluded  from  the  theory  because  the  heat  transfer  rates  are  zero  for 
high  blowing  rates . 

Additional  extensions  of  the  theory  might  entail  (1)  calculation 
of  the  inviscid  flow  past  the  point  where  the  injection  has  stopped, 

(2)  calculation  of  the  viscous  corrections  for  a  finite  Reynolds  number 
(which  will  require  a  pressure  correction  due  to  a  viscous  displacement 
thickness),  and  (3)  estimates  of  aerodynamic  coefficients  for  bodies  at 
small  angles  of  attack. 

An  unsolved  basic  problem  of  the  elementary  theory  is  whether  solu¬ 
tions  exist  for  constant  blowing  or  constant  pressure  or  both.  Pre¬ 
sumably  this  type  of  flow  in  two  dimensions  demands  a  more  accurate 
theory.  One  example  for  flow  over  a  cone  with  constant  blowing  is 
given  in  Ref,  11.  When  the  assumptions  of  the  present  theory  are  sat¬ 
isfied,  the  examples  discussed  in  the  preceding  sections  can  be  compared 
with  the  results  of  other,  more  elaborate  numerical  analyses.  Compari¬ 
son  of  theory  with  experiment  is  not  yet  possible,  since  the  injection 
distributions  of  the  present  work  have  not  yet  been  simulated  in  wind- 
tunnel  tests. 
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Append  ix 

SOLUTION  OF  THE  INTEGRAL  EQUATION  (2.30) 


The  basic  integral  equation  to  be  solved  in  the  region  x  ^  0  is 
(2.30): 


I/O  x  v  (5) 

(2a)1/2S(x)  =  N - - - 

0  (P(5)  -  P(x) 


(A-l) 


Equation  (A-l)  is  transformed  into  an  Abel- type  equation  by  regarding 
P(0)  -  P(x)  as  a  new  variable.  Let 


M.  =  P(0)  -  P(|) 


T1  =  P  (0)  -  P(x) 


(A-2) 


assuming  that  P(0)  is  finite  and  that  |j,,T|  are  monotone  functions.  Then 


(2a)1/2s(x(Tl))  =  } 


o  ('n  -  m-)1 


(A- 3) 


In  this  form  the  equation  is  a  standard  Abel  type  (Ref.  10)  and  has  the 
solution 


-Wi  ■ 


2o)1/Z  d_  j!  s(g(n) 

TT  dT|  j 


o  on  -  m-)' 


(A- 4) 


Integration  by  parts  and  the  use  of  S(0)  =  0  thus  yields 


,M  £  -  i 


o  on  -  n)' 


(A- 5) 
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The  reintroduction  of  the  original  variables  of  Eq.  (A-2)  gives 
the  final  result: 


v  (x)  = 

w'  ' 


_  i -2 22 


1/2 


TT 


©I 


S  (5) 


o  (p(5)  -  P(x) 


1/2 


d§ 


(A- 6) 


which  is  identical  to  Eq.  (2.31). 
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